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Abstract 

The Skyrme energy density functional has been applied to the study of heavy-ion fusion reactions. 
The barriers for fusion reactions are calculated by the Skyrme energy density functional with proton 
and neutron density distributions determined by using restricted density variational (RDV) method 
within the same energy density functional together with semi-classical approach known as the 
extended semi-classical Thomas-Fermi method. Based on the fusion barrier obtained, we propose 
a parametrization of the empirical barrier distribution to take into account the multi-dimensional 
character of real barrier and then apply it to calculate the fusion excitation functions in terms of 
barrier penetration concept. A large number of measured fusion excitation functions spanning the 
fusion barriers can be reproduced well. The competition between suppression and enhancement 
effects on sub-barrier fusion caused by neutron-shell-closure and excess neutron effects is studied. 
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I. INTRODUCTION 

A large number of fusion excitation functions have been accumulated in recent decades Q, 
j^, which provides a possibility for a systematic study on fusion reactions. Newton et.al. 
analyzed a total of 46 fusion excitation functions at energies above the average fusion 
barriers using the Woods-Saxon form for the nuclear potential in a barrier passing model of 
fusion. They found that the empirical diffuseness parameters a ranging between 0.75 and 1.5 
were considerably larger than those obtained from elastic scattering data and the deduced a 
showed strong increase with increasing the charge product Z\Z2- Thus, it results in a certain 
difficulty for giving satisfied predictions of fusion cross sections for unmeasured reaction 
systems. The fusion coupled channel model is successful for describing fusion excitation 
function of heavy-ion reaction at energies near fusion barrier. However, with the increasing 
of the neutron excess and the charges (the product Z\Z<i) of two nuclei, the fusion coupled 
channel model encounters a lot of difficulties due to a very large number of degrees of freedom 
involved. In order to carry out a systematic study of fusion excitation functions, a simple 
and useful approach seems to be required. A semi-empirical approach in which the quantum 
penetrability of Coulumb barrier is calculated by using the concept of barrier distribution 
arising due to the mult i- dimensional character of the real nucleus-nucleus interaction j^O, Q 
is very helpful in this aspect. Thus, a reasonable parametrization of the weighing function 
describing the barrier distribution based on the interaction potential in the entrance channel 
seems to be very useful. 

In (sl the interaction potential in the entrance channel for fusion reactions was calculated 
based on the semi-classical expressions of the Skyrme energy density functional js, 10, 11|. In 
the calculation of the interaction potential, which is the difference between the energy of the 
total system and the energies of individual projectile and target, the density distributions 



of th e p roj ee til e - t a ,et are reqU ,e, In Q the d e nSlty distri hu tionS of the pr0J ee t „e 
and target were determined by Hartree-Fock-Bogoliubov calculations. According to the 
Hohenberg and Kohn theorem[12], the energy of a N-body system of interaction fermions 
is a unique functional of local density. In the framework of the semi-classical Extended 
Thomas Fermi (ETF) approach together with a Skyrme effective nuclear interaction such 
a functional can be derived systematically. The density functional theory is widely used 
in the study of the nuclear ground state which provides us with a useful balance between 
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accuracy and computation cost allowing large systems with a simple self-consistent manner. 
Thus, it is very suitable to perform the calculation of the entrance-channel potential and 
the densities of reaction partners based on the same Skyrme energy density functional. Of 
course, the different density distributions obtained by different approaches will influence the 
entrance-channel potential. While a self-consistent treatment for both density distribution 
and the entrance-channel potential seems to be more reasonable. 

In this work, the Skyrme energy density functional is applied to make a systematic study 
of fusion reactions. Firstly, we will use the semi-classical expressions of the Skyrme energy 
density functional to study the energies and the density distributions of a series of nuclei 
by the restricted density variational (RDV) methodflO, [ill 13, 13 • Secondly, with the 
density distributions obtained, the entrance-channel potentials of a series of fusion reactions 
are calculated. Then, based on the entrance-channel potential obtained, a parametrization 
of the empirical barrier distribution is proposed to take into account the multi-dimensional 
character of real barrier and then apply it to calculate the fusion excitation functions in terms 
of barrier penetration concept. The paper is organized as follows: In Sec. II, the properties 
of ground state of nuclei and the entrance-channel potential are studied in the framework of 
the Skyrme energy density functional. In Sec. Ill, an approach to calculate fusion excitation 
functions is introduced and a large number of calculated results are presented. Finally, the 
summary and discussion are given in Sec. IV. 

II. THE CALCULATIONS OF FUSION BARRIER 

The energy density functional theory is widely used in many-body problems. In the frame- 
work of the semi-classical Extended Thomas Fermi (ETF) approach together with a Skyrme 
effective nuclear interaction, the energy density functional can be derived systematically. We 



take the Skyrme Hartree-Fock (Skyrme HF) formalism of the energy density functional[9, 
and based on it we calculate the proton and neutron densities of nuclei by means of restricted 
density variational methodjiol. 11, K^Ql- With the neutron and proton densities deter- 
mined in this way we calculate the fusion barrier for fusion reaction based on the same 
energy density functional. 



A. Skyrme Energy Density Functional 
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The total binding energy of a nucleus can be expressed as the integral of energy density 
functional 



E 



H dr. 



(1) 



The energy density functional Ti includes the kinetic, nuclear interaction and Coulomb 
interaction energy parts 

h 2 



n = ^M r ) + r n (r)] + H sky (r) + H coul (r). 



(2) 



For the kinetic energy part, the extended Thomas- Fermi (ETF) approach including all terms 
up to second order in the spatial derivatives (ETF2), is applied as that was done in ref.jsj]. 
With the effective- mass form factor 
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the kinetic energy densities r for protons (i = p) and neutrons (i = n) are given by 

5 36 pj 3 6 fi 
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where p^ denotes the proton or neutron density of the nucleus and p = p p + p n , Wo denotes 
the strength of the Skyrme spin-orbit interaction. The nuclear interaction part with Skyrme 
interaction H s k y reads 
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skyir) = ^ [(1 + ^x )p 2 - (x + -)(pp + p 2 )] 
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where to, ti, t 2 , £3, Xq, x±, x 2 , £3, and a are Skyrme-force parameters |8l. Ill | . The last term in 
the right hand of expression (5) is the semi-classical expansions (up to second order in h) of 
spin-orbit densities |9j. The Coulomb energy density can be written as the sum of the direct 
and exchange contribution, the latter being taken into account in the Slater approximation, 



(6) 



From the above expressions (l)-(6), one can see that the total energy of a nuclear system 
can be expressed as a functional of protons and neutrons densities [p p (r), p n { r )} under the 
Skyrme interaction associated with the ETF approximation. 

B. The Neutron and Proton Densities of Nuclei 

By minimizing the total energy of the system given by expression (1), the neutron and 
proton densities can be obtained, that is to solve the variational equation 
S f 

J {^K(r), p p (r)] - X n p n (r) - X P p P (r)} dr = 0, (7) 

with the Lagrange multipliers A n and X p to ensure the conservation of neutron and proton 
number. This density variational problem has been solved in two different ways in the past: 



13. 



14j resulting from Eq.(7) or by carrying 



±iQ. ( ( ) or by carrj 

m y Q Q- 



either by resolving the Euler-Lagrange equation 

out the variational calculation in a restricted subspace of functions HU Hl|, 112 I, 14 . In 
this work we take the neutron and proton density distributions of a nucleus as spherical 
symmetric Fermi functions 

-i 



pi(r) = p 0i 



1 + exp 



0/ 



{n,p} . 



(8) 



For the three quantities poi , Roi and cij in the equation, only two of them are independent 
because of the conservation of particle number = J pi(r)dr, iVj = {iV, Z}. For example, 
p 0p can be expressed as a function of R 0p and a p , 

-l 
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1+Tl 2 



R 



Op 



(9) 



with high accuracy when R 0p 3> a p . Here, R 0p , a p R 0n , the radius and diffuseness 

for proton and neutron density distributions, respectively. 

By using optimization algorithm, one can obtain the minimal energy Ey, as well as the 
corresponding R 0p , a p , R 0n , a n for the neutron and proton density distributions. 



The Skyrme force SkM*^?| is adopted in the calculations since SkM* is very successful 
for describing the bulk properties and surface properties of nuclei. It is a well known fact 
that ETF calculations with a reasonable effective interaction reproduce experimental bind- 
ing energies with a very good accuracy, when Strutinski-type shell corrections and pairing 
corrections are taken into account. And the obtained charge root-mean-square radii in this 
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work for shell closed nuclei are very close to the corresponding experimental data 
In addition we have made a comparison between the results with RDV method and the 
Skyrme-HF method and find that the surface diffuseness obtained by RDV method are little 
smaller than that by Skyrme-HF calculations. It is due to the neg lect of the higher order 
term corrections in the extended Thomas- Fermi (ETF) approach |10|. We have made a check 
that the surface diffuseness calculated with RDV method increases if all terms up to fourth 
order in the spatial derivatives in ETF energy density functional (ETF4) are taken into 
account. 

C. the Calculations of Fusion Barriers 

The interaction potential Vb(R) between reaction partners can be written as 

V b (R) = E tot (R)-E 1 -E 2 , (10) 

where R is the center-to-center distance between two nuclei, the E to t(R) is the total energy 
of the interacting nuclear system, E\ and E2 are the energies of individual nuclei (projectile 
and target), respectively. The interacti on p otential Vb(R) is also called entrance-channel 



on p c 



potential in ref.|a| or fusion barrier 2JJ, |22j, in the following we take the term of fusion 
barrier. The E to t{R), E±, E 2 are calculated with the same energy-density functional as that 
is used in the calculations of nuclear densities, 

Etot(R) = J U[pi p {r) + p 2p (r - R),pi„(r) + p 2n (r - R)] dr, 

Ei = J W[pi p (r),p ln (r)]dr, (11) 
E 2 = J H[p 2p (r),p 2n (r)]dr. (12) 

Here, pi p , p 2p , p\ n and p 2n are the frozen proton and neutron densities of the projectile and 
target, determined in the previous section. 




FIG. 1: The entrance-channel fusion barrier for 28 Si+ 92 Zr 
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For calculating the nuclear interaction energies and Coulomb energies, multi-dimensional 
integral 2^ is performed. For a certain reaction system, the fusion barrier is calculated in a 
region from R = 7fm to 20/m with step of AR = 0.25 fm. Fig.l shows the fusion barrier for 
the reaction 28 Si+ 92 Zr. The solid curve denotes the results based on the density distributions 
obtained by RDV method with ETF2. The dotted curve denotes the proximity potential 25}. 
The results with the density determined by restricted density variational calculations are 
very close to those of proximity potential at the region where the densities of two nuclei do 
not overlap. 

The fusion barriers for more than 80 fusion systems (ZiZ 2 > 150) are calculated based 
on ETF2. Part of them are listed in Tabled. Tabled lists the barrier height Bo, radius Rq 
(see Fig.l) and the curvature fruo of the barrier for a series of reactions. Here the curvature 
of the barrier, ftwo, is obtained approximately through fitting the barrier at the region from 
Rq — 1.25/m to Ro + 1.25/m by an inverted parabola. From the table, one can find that 
the curvature of the fusion barrier increases and the Q-value for complete fusion reaction 
decreases with the increase of product Z\Z 2 . 

For exploring the influence of order h 4 terms jlfli \v\ . full fourth order ETF (ETF4) 
is applied to calculate both the densities and the fusion barriers for reactions 28 Si+ 28 Si, 
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FIG. 2: The fusion barriers for 28 Si+ 28 Si, 40 Ca+ 48 Ca, 16 O+ 208 Pb and 48 Ca+ 208 Pb. The dotted 
curves and the solid curves denote the results with ETF2 and ETF4 approach, respectively, and 
the crossed curves and the dashed curves denotes the results of proximity potential and those of 
the analytical form proposed in 2^] , respectively. 

40 Ca+ 48 Ca^ 6 O+ 208 Pb, and 48 Ca+ 208 Pb. The similar work had been done by A. Dobrowolski 



et al. in 



lOt 111 | are included in the determination 



22]. We find that when the h A terms 
of the neutron and proton densities, the surface diffuseness of densities of nuclei increases 
and thus the height of fusion barrier calculated is lower (about 1 ~ 2MeV) than those with 
ETF2. The fusion barriers for these reactions calculated with different approach are shown 
in Fig.2. The dotted curves and the solid curves denote the results with ETF2 and ETF4 
respectively, and the crossed curves and the dashed curves denote the results of proximity 
potential and those of the analytical form proposed inj^], respectivel y. Our calculation 
results of barriers with ETF4 are very close to the analytical form int22]- We find that 
the barriers calculated with ETF2 are more close to the proximity potential for light and 
medium-heavy systems (for example ZiZ 2 < 680), while the barrier calculated with ETF4 
are more close to those of proximity potential for heavy systems. It implies that including 
fourth order terms of ETF approach is important in improving the surface diffuseness for 
heavy nuclei. 
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FIG. 3: a) The barrier heights -Bo f° r a series of reaction systems and b) The relative deviations 
from the results of proximity potential. The circles and crosses in Fig. 3(a) denote the calculation 
results with RDV and proximity potential, respectively. The solid line is the linearly fitting to our 
calculations. 

Fig. 3 shows the barrier heights Bq of the selected systems with ETF2 approach and the 
relative deviations of B Q from the proximity potentials i? prox [25] as function of ZiZ 2 . The 
circles and crosses in Fig. 3 (a) denote the results of our calculations and those of proximity 
potential, respectively. Our calculation results for the barriers are quite close to the prox- 
imity potential and the relative deviations are less than 2.5% (see Fig.3(b)) generally. From 
the figure one can also see that the barrier heights linearly increase with product Z\Z 2 - The 



solid line in Fig.3(a) is a linear function fit of the calculation results, i.e. Bju 
with R fit = 13.7 fm and V = 13.0MeV. 
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III. FUSION EXCITATION FUNCTIONS 



A. The Parametrization of Barrier Distribution 



According to Wong's formula 
barrier can be expressed as 



261 ] . the fusion excitation function for penetrating a parabolic 
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TABLE I: The fusion barriers of a series of reactions with ETF2 approach. 



Reaction R {fm) B (MeV) hw (MeV) Z 1 Z 2 Q-value(MeF) Reference 



28 Si+ 28 Si 


8.75 


29.86 


2.37 


196 


10.91 


[37] 


12 C+ 92 Zr 


9.75 


33.34 


2.10 


240 


0.94 


n] 


16 O+ 70 Ge 


9.25 


36.52 


2.84 


256 


2.51 


m 


16 0+ 72 Ge 


9.5 


36.27 


2.37 


256 


6.30 


m 


16 0+ 73 Ge 


9.5 


36.16 


2.44 


256 


8.84 


m 


16 0+ 74 Ge 


9.5 


36.08 


2.50 


256 


10.61 


m 


16 0+ 76 Ge 


9.5 


35.80 


2.63 


256 


10.51 


m 


12 C+ 128 Xe 


10.25 


40.98 


2.50 


312 


-0.91 


[39] 


16 + 92 Zr 


9.75 


43.88 


2.86 


320 


-3.94 


n] 


i6 0+ H2 Cd 


10.0 


51.12 


3.19 


384 


-9.91 


[40] 


27 Al+ 70 Ge 


9.75 


57.37 


3.26 


416 


-5.17 


[43] 


27 Al+ 72 Ge 


9.75 


57.05 


3.39 


416 


-4.21 


[43] 


27 Al+ 73 Ge 


9.75 


56.87 


3.46 


416 


-2.90 


[43J 


27 Al+ 74 Ge 


9.75 


56.72 


3.51 


416 


-3.21 


[43J 


27 Al+ 76 Ge 


10.0 


56.31 


3.01 


416 


-2.39 


[43] 


35 Cl+ 54 Fe 


9.5 


61.89 


3.67 


442 


-17.77 


[41] 


i6 0+ i44 Nd 


10.5 


61.35 


3.27 


480 


-22.43 


[45] 


12 C+ 204p b 


11.0 


59.80 


3.19 


492 


-28.40 


[52] 


16Q + 144g m 


10.5 


63.62 


3.27 


496 


-28.55 


12} 


17 0+ 144 Sm 


10.5 


63.26 


3.41 


496 


-24.43 


\2] 


16 0+ 147 Sm 


10.5 


63.31 


3.37 


496 


-24.65 


[47] 


16 0+ 148 Sm 


10.5 


63.20 


3.41 


496 


-23.09 


[2J 


16 0+ 149 Sm 


10.5 


63.08 


3.46 


496 


-21.71 


[47] 


16Q + 150g m 


10.5 


62.96 


3.49 


496 


-20.21 


[47] 


16 0+ 154 Sm 


10.75 


62.43 


3.07 


496 


-16.43 


\2] 


16 0+ 166 Er 


10.75 


67.83 


3.45 


544 


-25.13 


m 


28 Si+ 92 Zr 


10.0 


74.20 


3.97 


560 


-28.12 
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Reaction R (fm) B (MeV) hw (MeV) Z X Z 2 Q-value(MeF) Reference 



16q + 186 W 


11.0 


72.26 


3.52 


592 


-21.30 


1 

J 


32 S+ 89 Y 


10.25 


82.31 


3.76 


624 


-36.58 


[501 


33 S+ 90 Zr 


10.25 


84.12 


3.87 


640 


-39.76 


[511 


33 S+ 92 Zr 


10.25 


83.74 


3.97 


640 


-35.51 


[5JJ 


16 O+ 208p b 


11.25 


78.46 


3.54 


656 


-46.49 


[3] 


35 C1+ 92 Zr 


10.25 


88.61 


4.14 


680 


-39.37 


ft] 

J 


19 F+ 197 Au 


11.25 


85.21 


3.69 


711 


-35.92 


[52] 


16 + 232 xh 


11.5 


84.47 


3.68 


720 


-36.53 


[53, 541 


32 S+ 110p d 


10.5 


94.00 


4.13 


736 


-35.37 


[551 


36 S+ 110p d 


10.5 


92.85 


4.28 


736 


-38.01 


[55] 


19 F+ 208p b 


11.25 


87.44 


3.87 


738 


-50.07 




40 Ca+ 90 Zr 


10.25 


103.74 


4.31 


800 


-57.27 


[57] 


40 Ca+ 96 Zr 


10.5 


102.22 


4.28 


800 


-41.13 


[57] 


so Ti+ 90 Zr 


10.5 


111.16 


4.49 


880 


-64.73 


[58] 

HUM 


32 S+ 154 Sm 


11.0 


120.22 


4.47 


992 


-60.69 


M 

HUM 


28 Si+ 178 Hf 


11.25 


120.65 


4.41 


1008 


-64.77 


[61] 


29 Si+ 178 Hf 


11.25 


120.31 


4.42 


1008 


-65.70 


[61] 


30 Si+ i86 W 


11.5 


122.14 


4.24 


1036 


-70.22 


[521 


31p + 175 Lu 


11.25 


126.97 


4.46 


1065 


-70.45 


[61] 


28 Si+ 198p t 


11.5 


128.04 


4.38 


1092 


-78.75 


[621 


32 S+ 181 Ta 


11.25 


138.15 


4.56 


1168 


-80.58 


Ml 


32 S+ 182 W 


11.25 


140.02 


4.52 


1184 


-84.93 


\m 


132 Sn+ 64 Ni 


11.5 


162.75 


4.61 


1400 


-111.05 


[65] 


40 Ca+ 208p b 


11.75 


186.57 


4.62 


1640 


-136.69 


[66] 


48 Ca+ 208p b 


12.0 


183.17 


4.43 


1640 


-153.80 


m 



where E c . m . denotes the center-of-mass energy, B , R and Hul> are the barrier height, radius 
and curvature, respectively. This expression is based on the one-dimensional barrier pene- 
tration model. The one-dimensional barrier penetration model with empirically determined 
potential parameters is successful in describing the fusion excitation functions for light sys- 
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FIG. 4: The fusion barrier distribution for 16 0+ 92 Zr. The distribution is evaluated with AE c m , = 
1.8MeV. The solid dots and solid curve denote the experimental data and our calculation results, 
respectively. 

terns and heavy systems at energies above the barrier with some exceptions, but fails in 
describing sub-barrier fusion for heavy systems. It is found that for sub-barrier fusion of 
heavy systems, the measured fusion cross sections are of up to several orders of magnitude 
higher than the predictions of the model. 

Taking into account the multi-dimensional character []| of the realistic barrier we may 
introduce the barrier distribution in order to calculate the total fusion cross sections. Then, 
the fusion excitation function in terms of barrier penetration concept is given by 



where ai^ s (E c _ m _, B) is the fusion excitation function for the barrier B based on one- 
dimensional barrier penetration, D(B) is a weighing function to describe the barrier dis- 
tribution. Here we replace the barrier height Bq in Wong's formula Eq.(13) by B. The 
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FIG. 5: The weighing function for the fusion reaction of 16 0+ 186 W. The dashed and dot-dashed 
curves denote Gaussian distributions D\{B) and D2(B), respectively. The solid curve denotes the 
average value, i.e. D avT (B) = \D\{B) + D2(B)]/2. The vertical dotted lines denote the positions 
of Bq and B c , respectively. The arrow B x denotes the cross point between the curves D\{B) and 
D?,{B) at the left side of the peak of D avr (B). 



distribution function D(B) satisfies 



D(B)dB = 1. 



(15) 



D{B) is often taken to be continuous and symmetric distributions of rectangular or Gaussian 



shapes [2 7, 



28, 



29|. The Wong's formula is a special case when a delta distribution (D(B) 



8(B — B )) is taken. 

In this work, we propose a parametrization of the weighing function based on the fusion 
barrier obtained in the previous section. That is that we try to propose a macroscopic 
empirical barrier distribution rather than explicitly taking into account the coupling of the 
fusion motion to internal degrees of freedom as that is done in the fusion coupled channel 
modeljscj. Concerning the weighing functions, we first investigate the shape of barrier 
distributions extracted from experiments. We find in most cases the barrier distributions 
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are not symmetric. For example, in Fig.4 solid circles show the experimental fusion barrier 
distribution of 16 0+ 92 Zr (3^. One can see that the fusion barrier distribution of the reaction 
system is asymmetric, its left side is steeper than the right side. It indicates that only 
one Gaussian distribution is not good enough to describe the fusion barrier distribution. 
Motivated by the shape of the barrier distribution extracted from experiments, we consider 
the weighing function to be a superposition of two Gaussian functions. Two Gaussian 
distributions D\(B) and D 2 (B) are proposed as 



and 



with 



D X {B) 



v/7 



exp 



-7- 



(B-Bj 
(2^i) 2 



D 2 {B) 



2^pKW 2 



cxp 



(B 



B,) 2 



Wi 



\(Bo 



B c ), 



(16) 



(17) 



(18) 



w 2 



\(B 



(19) 



B x = B c + wi, 



(20) 



B 2 = B C + w 2 . 



(21) 



The barrier height B , radius Rq and curvature hhj Q are obtained from Tabled. The B c = 
fB is the effective barrier height taking into account the coupling effects to other degrees 
of freedom, such as dynamical deformation, nucleon transfer etc. We take the reducing 
factor / = 0.926. The 7 in -Di(-B) is a factor which influence the width of the distribution 
Di(B). The larger the 7 value is, the narrower the distribution is. For the fusion reactions 
between nuclei with non-closed-shell but near the /3-stability line, we take 7 = 1; for the 
reactions with neutron closed-shell nuclei or neutron-rich nuclei the value of 7 is calculated 
by a parameterized formula which will be discussed in the following sub-section. From the 
expression of (16) and (17) one can find that the peaks and the widths of D\(B) and D 2 (B) 
only depend on the height of the fusion barrier B except the 7 in Di(B). The peaks of two 
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Gaussian distributions locate at different energies between B c and B . For light systems, the 
width of the distribution becomes very small and thus it is very close to delta distribution, 
in consistent with the one-dimensional penetration model. Both Di(B) and D 2 (B) satisfy 
Eq.(15). Further we introduce D aYT (B) = (D 1 (B)+D 2 (B))/2. D avr (B) also satisfies Eq. (15). 
In Fig. 5 we show the parameterized weighing function of 16 0+ 186 W. The dashed and dot- 
dashed curves denote Gaussian distributions D\(B) and D 2 (B), respectively. The solid curve 
denotes the D mr (B). There are two crossing points between D 1 (B) and D 2 (B) (or D avI (B)) 
at both sides of D avI (B). The left one is located at B x (see Fig. 5). We notice that the left 
side of -D avr (-B) is too flat and does not fit the shape of the barrier distribution extracted 
from experimental data well (for example see the shape shown in Fig. 4). So we propose an 
effective weighing function 

f DAB) : B < B x 
D cS (B) = { (22) 

[D avI (B) : B>B X 

(with J D eS (B) dB m 1). We notice that the distribution Di(B) plays a main role for 
sub-barrier fusion and -D avr (.E>) contributes more to the fusion at energies near and above 
the barrier. 

With the D e ff(B) we can, in principle, calculate the fusion cross sections by the expression 

POO 

<Jfus(E c . m .)= / D cS (B)afl(E c . m .,B)dB. (23) 
Jo 

However, the D e g(B) does not satisfy Eq.(15) exactly because there is a small difference 
between D 1 (B) and D^^B) in the left side of the B x (see Fig. 5). To remedy this defect, we 
replace the expression (23) by the the following expression in the calculation of fusion cross 
sections 

Vf us {E c . m .) = min[CTi(£ c . m .),<7 avr (£ c . m .)]. (24) 

Here the <7i(.E c . m .) and o" avI .(-E c .m.) are calculated by expression (14) with D(B) = Di(B) 
and D(B) = D avr (B), respectively. Now the normalization condition Eq.(15) is always 
satisfied for the weighing function adopted in the calculations of fusion cross sections at 
each -E c . m .. We have checked that the fusion cross sections calculated with the expression 
(24) are very close to those calculated with the expression (23) when the 7 is not too 
large. In fact, we find that the results calculated with Eq.(24) is in better agreement with 
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FIG. 6: The fusion excitation functions for 16 0+ 186 W and 28 Si+ 92 Zr. The dashed and the dot- 
dashed curves denote the cross sections cri(£' c . m .) and cr avr (£' c . m .) calculated with D(B) = D\(B) 
and D(B) = D avr (5), respectively. 

experimental data. To check the barrier distribution obtained with our model, we also show 
the barrier distribution obtained from the calculation in Fig.4 by the solid curve. The good 
agreement between experimental data and calculation results indicates the parametrization 
of the weighing function is quite reasonable. 

B. Fusion Reactions between Nuclei with Non-closed-shell but near /5-Stability 

Line 

In this sub-section, the fusion excitation functions for about 30 reaction systems are inves- 
tigated. For these reaction systems, the projectile and target nuclei are non-closed-shell but 
near the ^-stability line (The shell-closure effects of 16 are neglected in this work). The 
shell effect and excess neutron effect are weak for this kind reaction systems. Fig. 6 shows 
the fusion excitation functions of 16 0+ 186 W and 28 Si+ 92 Zr. From the figure one can find 
that at energies above the barrier a avr (dot-dashed curves) is more close to the experimental 
data, while at energies below the barrier, <Ji (dashed curves) is preferable. From Fig.6(c) 
and Fig. 6(f) one can see that all experimental data spanning the fusion barrier for both 
16 0+ 186 W and 28 Si+ 92 Zr can be reproduced well. Fig.7 to Fig.9 show more calculation re- 
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E c. m .( MeV ) E c.™.( MeV ) E c. m .( MeV ) E a m .( MeV ) 

FIG. 7: The fusion excitation functions for 12 C, 16 0, 33 S, 35 Cl+ 92 Zr and 16 0+ 112 Cd, 144 Nd, 166 Er, 
232 Th. The squares denote the experimental data and the solid curves denote the calculation results 
with 7=1. 

suits and experimental data for comparison for this kind fusion reactions. All the fusion 
excitation functions can be reproduced very well, which indicates our parametrization of 
Di(B) and D 2 (B) is quite useful and reasonable. 



17 



1000 
100 



E 



0.1 

0.01 



1000 
100 




(a) 

75 90 105 120 

E (MeV) 




100 1 20 140 1 60 1S0 




20 30 40 50 

E (MeV) 



1000 
100 



f 10 



a Si + ™Hf, 






■ exp. 




I calc. 




,(e) 



E c . m .(MeV) 




- ^S-^Ta 








■ exp. 
cdc. 




(h) 



110 1 20 130 1 40 150 

E cm (MeV) 



1000 
100 



10 




100 110 1 20 1 30 140 150 

E (MeV) 



1000 
100 



10 



30 Si+ 136 W_ 






■ ex[j. 




calc. 




(f) 



100 110 120 130 140 130 



105 120 135 150 165 

E (MeV) 



1000 
100 



E 



, 32 c 182,., _^ 








■ exp. 




calc. 




(i) 



120 140 160 180 200 

E (MeV) 



120 140 160 180 200 

E cm (MeV) 



FIG. 8: The fusion excitation functions for 19 F+ 197 Au, 28 Si+ 28 Si, 28 ' 29 Si+ 178 Hf, 28 Si+ 198 Pt, 
28 Si+ 186 W, 31 P+ 175 Lu, 32 S+ 181 Ta and 32 S+ 182 W. The squares denote the experimental data and 
the solid curves denote the calculation results with 7=1. 







16 0+ 73 Ge^_— 


i- Xaj^. 




/ I ' exp \ 

1 1 calc.1 


/■ I ■ exp. 1 
r I c^Jc.l 


/ I ■ exp. 1 
/ I calcl 


/ I ■ exp. 
-, - J I calc. 


/ I ■ exp. I 
/ I calcl 


,/ (a) 


7 (b) 


,/ (c) 


,/ W) 


/ , , , , m 


30 35 40 45 50 


30 35 40 45 50 i 


530 35 40 45 50 


5530 35 40 45 50 55 


30 35 40 45 50 55 


E c. m .< M eV) 


E c. m .< M eV) 


E „ m .< M eV) 


E c. m . 0*v) 


E , m .(MeV) 


^AI+^Ge^^- 


^AI+^Ge 


^AI-^Ge 






/ 1 ■ exp.| 
r 1 calcl 


/ I ■ exp. 1 
/■ I calc.li 


/ I ■ exp- 1 

/ I calcl 


/ I ■ exp-l 
A I calcl 


/ I ■ exp. 1 

/ I calcl" 


/ . (0 


/ (9) 


(h) 


/ . (0 


/ (i) 


48 52 56 60 64 48 52 56 60 64 48 52 56 60 6 


4 48 52 56 60 6 


X 48 52 56 60 64 


E .,,( M eV) 


E c.n,. < MeV > 


E.,„.(MeV) 


E c . m . (MeV) 


E cm . (MeV) 



FIG. 9: The fusion excitation functions for 16 O, 27 Al+ 70 ' 72 ' 73 ' 74 ' 76 Ge. The squares denote the 
experimental data and the solid curves denote the calculation results with 7 = 1. 
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TABLE II: The AQ and 7 values for a series reactions with neutron closed-shell nuclei or neutron- 
rich nuclei. 



Reaction 


AQ(MeV) 


7 


Reference 


33 S+ 90 Zr 


-4.25 


3.6 


[5JJ 


50 Xi+ 90 Zr 


-3.25 


3.6 


[58] 


35 Cl+ 54 Fe 


-6.58 


4.8 


[41] 


16 O+ 208p b 





1.5 


13} 


19 F+ 208p b 





1.5 


m 


17 0+ 144 Sm 


-0.21 


1.6 


[2] 


32 S+ 154 Sm 


9.80 


0.5 


[59] 


32 S+ no pd 


4.13 


0.6 


[55] 


40 Ca+ 96 Zr 


10.00 


0.5 




132 Sn+ 64 Ni 


7.50 


0.7 


M\ 


16 0+ 144 Sm 


-3.91 


3.5 


12] 


16 0+ 147 Sm 





1.0 


[47] 


16 0+ 148 Sm 


1.56 


0.8 


[2] 


16 0+ 149 Sm 


2.94 


0.7 


[47J 


16 O+ 150g m 


4.44 


0.6 


[47J 


16Q + 154g m 


8.22 


0.5 


[2] 


16 0+ 144 Sm 


-3.91 


3.5 


[2] 


16Q + 154g m 


8.22 


0.5 


\2] 


28 Si+ 58 Ni 


-6.81 


4.4 


\m 


28 Si+ 64 Ni 


8.56 


0.5 


[67] 


36 S+ 90 Zr 


-1.25 


2.6 


[68] 


36 g + 96 Zr 


7.85 


0.7 


m 


40 Ca+ 90 Zr 


-6.14 


5.1 




40 Ca+ 96 Zr 


10.00 


0.5 


[57] 


32 S+ no pd 


4.13 


0.6 




36 S+ no pd 


1.49 


1.3 




40 Ca+ 48 Ca 


10.64 


0.9 


[42] 


48 Ca+ 48 Ca 


3.09 


1.7 


[42] 
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C. Fusion Reactions between Nuclei with Neutron-Shell-Closure or 

Neutron-rich Nuclei 



The shell effect for neutron closed-shell nuclei and the excess neutron effect for neutron-rich 
nuclei play an important role in the fusion reactions at energies below the fusion barrier. 
Fig. 10(a) shows the fusion excitation function of 33 S+ 90 Zr. The dashed curve denotes the 
calculation results with 7 = 1 and the squares denote the experimental data. One can find 
that at the energies near and above the barrier B the calculation results with 7 = 1 are 
in good agreement with the experimental data. However, at sub-barrier energies the fusion 
cross sections are over-predicted with 7 = 1. Through adjusting the 7 value, we find that the 
fusion excitation function at sub-barrier energies for 33 S+ 90 Zr can be reproduced reasonably 
well when 7 = 3.6. Fig.lO(b) shows the fusion excitation function of the neutron-rich 
nuclear fusion reaction 32 S+ 154 Sm. For this reaction, the fusion cross sections at sub-barrier 
energies are under-predicted when 7 = 1. For this case, a small 7 should be taken so that the 
distribution -Di(-B) is broaden and thus the fusion cross sections are enhanced at sub-barrier 
energies. We find that when 7 = 0.5 the sub-barrier fusion cross sections can be reproduced 
well (see Fig.lO(b)). Through systematically analyzing the fusion excitation functions of 
reactions between nuclei with neutron closed-shell but near the /3-stability line and those 
with neutron-rich nuclei, we find that fusion cross sections are suppressed for the former 
cases and enhanced for the latter cases compared to the calculation results with 7 = 1. The 
enhancement of fusion cross sections for reactions with neutron-rich nuclei compared with 
non- neutron-rich nuclei has been found in refs. ?| 31, 32 ]. It is attributed to the neutron 
transfer and neck formation which lower the fusion barrier and thus enhance the fusion cross 
sections at sub-barrier energies}?!.!^. For neutron closed-shell nuclei, the strong shell effect 
suppresses the lowering barrier effect. Based on above discussion we propose an empirical 
formula for the 7 values used in the weighing function -Di(-B) for systems with the same Zi 
and Z 2 , 

1 = l-c AQ + 0.5(5r 9 + 5 t D, (25) 

where AQ = Q — Q denotes the difference between the Q-value of the system under 
considering for complete fusion and that of the reference system. The Qo is the Q-value 
of reference system. Co = 0.5MeV~ l for AQ < cases and Cq = O.lMeV" 1 for AQ > 
cases. 8% , °i( ta rg) _ ^ £ Qr neu ^ ron closed-shell projectile (target) nucleus and ^°^ tar9 ^ = for 
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E, m .(MeV) E c . m .( MeV ) 
FIG. 10: The fusion excitation functions for reactions 33 S+ 90 Zr, 32 S+ 154 Sm. The squares denote 
the experimental data. The solid and dashed curves denote the calculation results with the 7 values 
obtained by formula (25) and with 7 = 1, respectively. 

non-closed cases (The shell-closure effects of 16 are neglected in this work as mentioned in 
the above sub-section). In addition, we introduce a truncation for 7 value, i.e. 7 should not 
be smaller than 0.5. 

The reference system is chosen to be the reaction system with nuclei along the /3-stability 
line. More precisely, we do it as follows: from the periodic table we find the relative atomic 
masses of corresponding elements of projectile and target (M^ r ^ , M* a ^ 9 ), then the mass 
numbers for reaction partners of reference system (A^° 3 and A 1 ^ 9 ) can be obtained by 
inequality A l — 1 < M* m < A % . % denotes projectile [i = proj) or target {i = targ) 
nuclei. For example, for three fusion reactions 16 O+ 90 ' 92 ' 96 Zr, 16 0+ 92 Zr reaction system 
is taken to be the reference system according to the atomic mass of O and Zr obtained 
from periodic table and inequality given above. If there exist certain reaction system for 
which the experimental data of fusion excitation function can be reproduced well by the 
calculations with 7 = 1, then this reaction system is preferably chosen as the reference 
system for the series of reactions with the same Z\ and Z 2 . For example, for the series 
of reactions of 16 0+Sm, the reaction 16 0+ 147 Sm is taken as the reference system because 
the experimental data can be reproduced by the calculations with 7 = 1 (see the following 
discussion). For this kind of reactions, there may exist possibility of more than one such 
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E =.m. ( MeV > E cm (MeV) E m (MeV) 

FIG. 11: The fusion excitation functions for reactions 33 S, 50 Ti+ 90 Zr, 35 Cl+ 54 Fe, 16 O, 19 F+ 208 Pb 
and 17 0+ 144 Sm with neutron-shell-closure nuclei near the /3-stability line. The squares denote the 
experimental data. The solid denote the calculation results with the 7 values obtained by formula 
(25). 

reaction systems for which the experimental data can be described by the calculations with 
7 = 1. This situation might be rare and from our investigation we have not encountered. If 
this situation occurs we prefer to choosing the reference system according to the inequality 
given above. The 7 values for some fusion reactions with neutron closed-shell nuclei or 
neutron-rich nuclei are calculated by expression (25) and the results are listed in Table. II. 
For fusion reactions between nuclei with neutron closed-shell and near the /^-stability line, 
the 7 values are larger than 1 (see Table. II), which means that the width of the barrier 
distribution D\(B) becomes narrow and thus the fusion cross sections of these systems at 
sub-barrier energies are suppressed compared with the reference system. Fig. 11 shows the 
fusion excitation functions of these fusion reactions. One can see from the figure that the 
experimental data are nicely reproduced. 

Further, we find that at energies near and above the fusion barrier, calculated fusion 
cross sections are not sensitive to the value of 7 (see Fig. 10 and Fig. 12(a) and Fig. 12(f)), 
which implies that we can calculate the fusion cross sections at energies near and above the 
fusion barrier for unmeasured fusion system by simply taking 7 = 1 in our parametrization 
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E m (MeV) E cm (MeV) E cm (MeV) 



FIG. 12: The fusion excitation functions for 16 Q+i44,i47,i48,i49,i50,i54g m S q Uares denote the 
experimental data. The solid and dashed curves denote the calculation results with 7 obtained by 
formula (25) and with 7 = 1, respectively. 

of weighing function. For sub-barrier fusion, more sophisticated investigation for 7 value is 
required. For further investigating the influence of 7 value, a series of fusion reactions with 
16 projectile on Sm isotope targets from 144 Sm to 154 Sm are studied. Fig. 12 shows the fusion 
excitation functions for 16 + 144_154 Sm. In the calculations we find that the fusion cross 
sections for 16 0+ 147 Sm can be reproduced well by the calculations with 7 = 1 (see Fig. 12(b)), 
so this system is taken as a reference system for studying other 16 0+ 144_154 Sm. The AQ 
for 16 0+ 144 Sm is negative as obtained from Table. II and consequently the corresponding 7 
value calculated by expression (25) is larger than 1. Thus, this fusion process is unfavored 
and the fusion cross sections at sub-barrier energies are suppressed compared with reference 
system 16 0+ 147 Sm. With the increase in numbers of neutrons, the Q- values of reaction under 
consideration increase gradually and the corresponding 7 values decrease, which indicates 
that the sub-barrier fusion cross sections change from suppression to enhancement compared 
with the reference system. Our calculated results for O+Sm reactions are in good agreement 
with the experimental data. For investigating the competition between suppression and 
enhancement effect of fusion cross sections, six pair reactions of 16 0+ 144 ' 154 Sm, 28 Si+ 58 ' 64 Ni, 
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FIG. 13: The fusion excitation functions for Sn+ Ni. The squares denote the experimental 
data and the solid curves denote the calculated results with 7 obtained by formula (25). The 



dot-dashed and doted curves denote the results of fusion coupled channel model 



351 ] . respectively. 



36g + 90,96 Zr; 40 Ca+ 90,96 Zr; 32,36g + 110p d and 40,48 Ca+ 48 Ca are studied . We find that those 

reactions with 7 < 1 are all enhanced at sub-barrier energies compared with non-neutron- rich 
systems, while those reactions with 7 > 1 are all suppressed at sub-barrier fusion compared 
with neutron non-closed-shell systems. 

For radioactive beam fusion reactions 132 Sn+ 64 Ni, the competition exists since 132 Sn 
is both neutron-rich and neutron-shell-closure (N = 82, Z = 50). The 7 value for this 
reaction is smaller than 1 (see Table. II) and the enhancement of the fusion cross sections 
at sub-barrier energies is expected. The fusion cross sections of this system are shown 
in Fig. 13. Squares and solid curve denote the experimental dataj^ and our calculation 
results, respectively. The results of fusion coupled-channel model (The dotted curve denotes 
inelastic excitations and the dashed curve denotes both the inelastic excitations and neutron 
transfer are considered 0]) are also presented for comparison. From the comparison one can 
find that the agreement in sub-barrier fusion cross sections calculated with our approach is 
better than the fusion coupled channel model [30] calculations. This seems to indicate that 
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there are still some physical aspects missing in present coupled channel calculations. 

IV. CONCLUSION AND DISCUSSION 

In this work, the Skyrme energy density functional has been applied to study heavy-ion 
fusion reactions. The properties of ground state nuclei are studied by using the restricted 
density variational method with the Skyrme energy density functional plus the semi-classical 
extended Thomas- Fermi approach (up to second order in K). With the proton and neutron 
density distributions obtained in this way, the fusion barriers of a series of reaction systems 
are calculated by the same Skyrme energy density functional. We propose a parametriza- 
tion for the weighing function describing the empirical barrier distribution based on the 
fusion barrier calculated with Skyrme energy density functional. The weighing functions 
of the barrier are assumed to be the superposition of two Gaussian functions. With the 
parametrization of the weighing function for the empirical barrier distribution, fusion exci- 
tation functions for more than 50 systems are calculated. A large number of measured fusion 
excitation functions spanning the fusion barriers can be reproduced well. The competition 
between suppression and enhancement effects on sub-barrier fusion caused by neutron-shell- 
closure and excess neutron effects have been investigated. 

However, the proton-sheU-ciosore effects, d_, effect, HQ as we., as the effects 
due to large deformation of nuclei in the fusion reactions have not been taken into account 
in present work, yet. All those effects are very important in fusion dynamics but they are 
beyond the scope of this work. The study on these aspects is under way. 
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